Thermopower of an SU(4) Kondo resonance under an SU(2) symmetry-breaking field 
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We calculate the thermopower of a quantum dot described by two doublets hybridized with two 
degenerate bands of two conducting leads, conserving orbital (band) and spin quantum numbers, as 
a function of the temperature T and a splitting 8 of the quantum dot levels which breaks the SU(4) 
symmetry. The splitting can be regarded as a Zeeman (spin) or valley (orbital) splitting. We use 
the non-crossing approximation (NCA), the slave bosons in the mean-field approximation (SBMFA) 
and also the numerical renormalization group (NRG) for large S. The model describes transport 
through clean C nanotubes and in Si fin-type field effect transistors, under an applied magnetic 
field. The thermopower as a function of temperature S(T) displays two dips that correspond to 
the energy scales given by the Kondo temperature Tk and 5 and one peak when ksT reaches the 
charge-transfer energy. These features are much more pronounced than the corresponding ones 
in the conductance, indicating that the thermopower is a more sensitive probe of the electronic 
structure at intermediate or high energies. At low temperatures (T <C Tk) TkS(T)/T is a constant 
that increases strongly near the degeneracy point 5 = 0. We find that the SBMFA fails to provide an 
accurate description of the thermopower for large S. Instead, a combination of Fermi liquid relations 
with the quantum-dot occupations calculated within the NCA gives reliable results for T <C Tk- 

PACS numbers: 72. 20. Pa, 72.15.Qm,73.23.Hk, 73.63.Nm 



I. INTRODUCTION 

Materials with potentially useful thermoelectric prop- 
erties and in particular large thermopower, are currently 
a subject of intense research due to their potential ap- 
plications, for example in refrigeration or conversion 
of waste heat into electricity^- In addition, the ther- 
mopower is a useful tool to obtain additional insight on 
fundamental problems, like the Kondo effect^— The fo- 
cus of research on thermal and thermoelectric properties 
has moved in the last years to nanostructures such as 
quantum dots?ii~— carbon nanotubes ) 17 ' 18 molecules^ 
nanowires, and spin systems^ - — As we will show, the 
thermopower can be a more useful tool than the conduc- 
tance, to study some features of the electronic structure 
of Kondo systems at finite energies. 

Electronic and thermal transport through single level 
quantum dots is well described by the SU(2) Anderson 
model, and the Seebeck coefficient (thermopower) has 
been calculated using this modeh 11 ^ 4 " — In particular, 
Costi and Zlatic have made a comprehensive study of 
the transport properties using the numerical renormal- 
ization group (NRG) A 1 - A particularly interesting mul- 
tilevel system is the SU(4) Anderson model, which de- 
scribes quantum dots in carbon nanotubes^^— and sili- 
con nanowires^ In contrast to the SU(2) case, the elec- 
tron spectral density near the Fermi level of the SU(4) 
Anderson model in the strong coupling (Kondo) limit for 
nearly one electron at the dot is highly asymmetric 1 24 i 28 
leading to a high derivative at the Fermi energy, which 
(neglecting phonon-drag effects) in turn is proportional 
to the Seebeck coefficient S at low temperatures 



The interest in Si nanowires increased due to the pos- 
sibility of reducing the thermal conductance k, lead- 
ing to a large figure of merit aS 2 T '/ k, where a is the 
conductivity^ Recently, the conductance of Si fin-type 
field effect transistors has been measured under an ap- 
plied magnetic field B, which leads to a crossover from 
an SU(4) to an orbital SU(2) Kondo effect 3 The results 
were interpreted using the non-crossing approximation 
(NCA). In C nanotubes, even for B — 0, there might be 
a symmetry reduction to SU(2) due to disorder-induced 
intervalley mixing^ 

In this work, we calculate the thermopower of the An- 
derson impurity model for two doublets, each one either 
spin or orbital degenerate, and infinite on-site Coulomb 
repulsion in the Kondo limit, as a function of the level 
splitting (5, which corresponds to the Zeeman splitting for 
two orbitally degenerate levels in presence of an applied 
magnetic field, as in Si fin-type field effect transistors^ 
The symmetry of the model is SU(4) for 8 — and SU (2) 
for 5^0. We use complementary theoretical approaches: 
the NCA^i^i 2 ^ slave bosons in the mean-field approxima- 
tion (SBMFA) j5Sr^& and Fermi liquid theory for low tem- 
peratures. We also use NRG to test the other approaches 
in the SU(2) limit. 

The SBMFA satisfies Fermi liquid relationships at 
T = and is expected to capture the low-energy physics 
in the Kondo limit. In turn, while NCA has an error of 
about 15% in the value of the spectral density at very low 
temperature according to the Friedel sum rule;* 8 , it de- 
scribes better the whole behavior, as shown for example 
in a previous comparison of NCA and NRG results in the 
SU(2) casei2£t In addition, if the differential conductance 
G is normalized at T = 0, the leading behavior of G for 
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small voltage and temperature 3 -! agrees with alternative 
Fermi liquid approaches ^Sr— and the temperature de- 
pendence of the conductance practically coincides with 
the NRG result over several decades of temperature^ 
Calculations of the thcrmopower for more than one spin- 
degenerate level within the NCA, compared well with 
experiments in some Ce compounds. 8 - 9 NRG is a very 
accurate technique at low temperatures.™ However, for 
two bands, the Hilbert space is increased 16 times in 
each iteration, instead of 4 times for one band, making 
the technique much more demanding if the same degree 
of accuracy is wished. In addition, the SU(4) symmetry 
cannot be used to reduce the size of the matrices, since 
this symmetry is broken by 5. Therefore we use here 
NRG only in the limit 6 — > +oo, in which only one dou- 
blet and the band that mixes with it play a significant 
role. 

The paper is organized as follows. In Section |TT] we 
describe the model used. The approximations and the 
equation for the Seebeck coefficient are described in Sec- 
tion IIIII The numerical results are presented in Section 
TVl Section|V]contains a summary and a short discussion. 



II. MODEL 

We start with a generalization of the Anderson model 
for infinite on-site Coulomb repulsion, which contains a 
singlet \s) with J\f (even) particles and two spin doublets 
icr) (i = 1, 2 is the valley index; a =t or \.) with J\f + 1 
(or A/" — 1) particles representing the four spin and valley 
degenerate states of a quantum dot created for example 
by depleting the density at two points of a C nanotube or 
in a Si nanowire. The dot is connected to two conducting 
leads which are also spin and valley degenerate. The 
SU(4) symmetry is then broken by applying a magnetic 
field or breaking the valley degeneracy in a simple way. In 
both cases, interchanging spin (cr) and valley (i) indices 
if necessary, the Hamiltonian can be written in the form 



H = E s \s)(s\+2_,Ei\ia)(i(T\ + 



vka 



We note that in the case of intervalley mixing in- 
duced by disorder in C nanotubes, the hopping elements 
for small magnetic field depend on the valley and lead 
indices ^2, In this case, the formalism used in this work is 
not applicable. It seems that a full non-equilibrium for- 
malism is needed to treat the most general that 
developed for the conductance in Ref. |43| and sketched 
in Ref. |3 The model is not applicable either for the 
case of magnetic impurities in C nanotubes, for which 
symmetry-breaking geometrical effects play an important 
role4£ 



III. THE FORMALISM 

A. Equations for the thermopower 

In the limit of vanishing applied bias voltage and tem- 
perature difference between the leads, the electronic part 
of the transport coefficients can be evaluated in terms of 
the total spectral density at the dot Pd{^>) = X) 

icr 

This is possible due to the fact that the couplings be- 
tween the quantum dot and the right or left leads are 
proportional. 46 If this were not the case (as in nanotubes 
affected by disorder— or for some molecules^ 3 -) a different 
formalism would be needed j 43 1 44 

The Seebeck coefficient is simply given by-ii 



S 



-h{T) 



(2) 



eTI (T)' 

where e is the absolute value of the electronic charge and 



/„ = / u) n p d (uj) ( ) du>, 



(3) 



where f(to) is the Fermi function, and the zero of energy 
is taken at the Fermi energy tp = 



+ ^{V v \ia){s\c vki<J +H.c), 



(1) 



where c vki create conduction states at the left (y = L) 
or right {v = R) lead, and V v is the hopping between 
the lead v and both doublets, assumed independent of 
k for simplicity. The doublets are split by an energy 
6 = E-2 — E\. This corresponds to the Zeeman splitting 
when the SU(4) symmetry of the model for S = is 
broken by an applied magnetic field. 

While there are four spin degenerate bands of mobile 
electrons, depending on valley index i or position with re- 
spect to the quantum dot (left or right), for each energy 
(Lk = £u,k' for which there are states at the left and the 
right, only the linear combination V L kc\ kilJ + V Rk >c Rk , ia 
hybridizes with the state \icr). Thus, the model is effec- 
tively a two-band model. 



B. Approximations for the spectral density 

To calculate the total spectral density Pd(w) that enters 
Eqs. © we have used mainly the NCA. At T = have 
used the SBMFA and also a combination of Fermi liquid 
relationships and quantum-dot occupations obtained us- 
ing NCA. In the limit 8 — > +oo, we have also used NRG 
to shed light on the virtues and shortcomings of the other 
approximations. 

In the NCA and SBMFA, an auxiliary boson b, and 
four auxiliary fermions fi„ are introduced, so that the 
localized states are represented as 



(4) 
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where |0) is the vacuum, 
satisfy the constraint 



These pseudoparticles should 



(5) 



The NCA solves a system of self-consistent equations 
to obtain the Green functions of the auxiliary particles, 
which is equivalent to sum an infinite series of diagrams 
(all those without crossings) in the corresponding per- 
turbation series in the hopping, and afterwards a projec- 
tion on the physical subspace of the constraint is made. 
The formalism of the NCA for this problem or similar 
ones is described in previous papers* ^ 24 ' 43 In particu- 
lar, the more general case of complex hoppings is treated 
in Ref. H^. Therefore, we do not give more details here. 
The application of NRG in the case of only one doublet 
(5 — > +oo) has also been explained before^ 

In the SBMFA, the boson operators are replaced by 
a number bo = (b), and the energy is minimized with 
respect to bo and a Lagrange multiplier A that enforces 
Eq. ([3]). Assuming a constant density of unperturbed 
conduction electrons p extending from — D to D and filled 
to the Fermi level cf = 0, a simple generalization of the 
case of one doublet) 33 ! 34 leads to the following change of 
energy after introduction of the impurity 



nia = (\ia){ia\) = — arctan ( — | , 



&8A/7T 



(w - E^ 2 + A 2 



(9) 



C. Fermi liquid theory 



For an SU(N) model which channel index j = 1 to 
N, and a simple symmetry breaking perturbation (like a 
generalized magnetic field) that does not mix channels, 
so that spin and channel are conserved, the Friedel sum 
rule relates the spectral density at the Fermi level e f with 
the number of displaced electrons for each channel^! The 
latter coincides with the occupation for a constant unper- 
turbed density of conduction electrons with a wide band 
O > A as we assume, where A (of the order of the 
Kondo temperature Tk) is the resonant level width of 
the quasiparticles [as in the SBMFA, Eq. 

Then one has^il 
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Ef 



D 2 



2Ei arctan 



Ei 



+ A(6g-1), 



(6) 



where Ei = E t + A, A = 6gA, and A = np(V£ + Vg) 
is the total half resonant level width (adding the contri- 
butions from both leads). Above D 3> A was assumed. 
Minimizing E with respect to A one obtains an equation 
that allows to relate E\ (or A) with half the quasipar- 
ticle level width A (which is of the order of the Kondo 
temperature) . After some algebra we obtain 



^ + (1 + /3 2 )A 2 



-1 1/2 



/?A-~ 
2' 



E x = 
/T 1 = tan 

Minimizing E with respect to A one obtains 




(7) 



1 (^ 1 2 + A 2 )((^ 1 + J) 2 + A 2 ) 

7T 11 D A 



Ei — Ei 



= 0, (8) 



and replacing E\ from Eq. ((TJ) in Eq. ((5|) an equation 
for the single unknown A is obtained, which we solve 
numerically. 

The occupation and the spectral density near the Fermi 
energy for each doublet are 



Pj^f) 



7tA 



sin 2 (7rnj). 



(10) 



In addition, the derivatives at cp are also known (using 
for example renormalized perturbation theory) 41 i 48 i 49 , 



_ g _| WP = n(2 mv )- s -. 



(11) 



From Eqs. ©, it is apparent that the SBMFA (in 
which we have chosen ep = 0) satisfies these relation- 
ships. 

Using Eqs. (flQ)) and (TTT1) and a Sommerfeld expansion 
in Eqs ©, the Seebeck coefficient Eq. © for T -> 
becomes 



S=- 



2-K 2 k B T J2ia sin 3 (7rn i(T ) cos(7rn icr ) 



3A 



(12) 



IV. NUMERICAL RESULTS 

Without loss of generality, we take ep = E s = 0, where 
tF is the Fermi level of the leads. For the numerical 
calculations, we assume a constant density of states per 
spin of the leads p — 1/(2D) between — D and D. We 
take the unit of energy as the total level width of both 
doublets r = 2 A = 2%p(Vl + Vg). The energy of both 
doublets is denoted as Ei — Ed, E 2 — Ed + 6. 
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FIG. 1. (Color online) Thermopower as a function of tem- 
perature for Ed = —4T and 8 — »■ +00. Squares: NRG, dots: 
NCA, dashed line: Lorentzian spectral density (see text). 



A. The limit of one doublet 

For 5 +00, the model is equivalent to the limit of in- 
finite on-site repulsion U of the simplest Anderson model, 
studied in detail before using the NRGi^i Here we com- 
pare results obtained with NRG with those of NCA and 
SBMEA to see the limitations of the different methods, 
which will be useful for the analysis of the general case. 

In Fig. [TJwe display the thermopower as a function of 
temperature in the Kondo regime. There is a dip near the 
Kondo temperature Tk due to the fact that the Kondo 
peak in the spectral density is slightly above the Fermi 
level ep, and a peak near the charge-transfer energy ep — 
Ed due to the corresponding peak below ep. We mean 
"near" as an estimation of the order of magnitude. For 
example, the maximum of the curve within NCA is at 
T = l.lOr, while e F - E d = AT in the figure. 

The spectral density pd(u>) of the Anderson model for 
only one doublet (S — > +00) is well known. In particu- 
lar, using a local-moment approach, Logan and coworkers 
have shown that in the Kondo regime, the charge-transfer 
peaks (which carry most of the spectral weight in this 
regime) have a Lorentzian shape of a width two times 
larger than that of the non-interacting case^ We have 
included in Fig. [T]the thermopower that results replac- 
ing Pd{u) in Eqs. (|2|) and (|3J|, by a simple Lorentzian of 
total width 2r centered at Ed- The agreement with the 
NCA result at high temperatures is remarkable. Instead, 
the charge-transfer peak in the thermopower within NRG 
has a small shift to higher temperatures and an intensity 
about 10 % lower. We have verified that this also hap- 
pens for other values of Ed- This is probably related 
with resolution problems of the NRG at large energies, 
as discussed below. 

At low temperatures, the NRG results are more reli- 
able than those of NCA and SBMFA. We have extracted 



TABLE I. Kondo temperature from G(T) (Tk), total occupa- 
tion of the lowest doublet (n), linear coefficient of S(T) from 
NRG (STg/T) and from a Fermi liquid theory (SA/T) for 
different values of the charge-transfer energy — Ed- 
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-3.48 
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2.22 x 10" 5 


0.941 
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0.958 


-0.38 


-0.437 
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3.16 x 1CT 8 


0.967 


-0.27 


-0.339 
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1.18 x 1(T 9 


0.973 


-0.22 


-0.317 



the linear term in T in the Seebeck coefficient, plotting 
S/T vs T for small T and looking at the extrapolation to 
T = 0. The results are shown in Table I. They have an 
uncertainty of the order of 10 %. We have also calculated 
the conductance G(T) (not shown) and the total occu- 
pancy n adding both spins. Using the results for G(T), 
we estimated the Kondo temperature from the require- 
ment that G(T%) = G(0)/2. Using n n = n n = n/2, 
712^ = ?i2j. = 0, and the Fermi liquid expression Eq. (|12p . 
S/T is calculated in an independent way. Taking into 
account the exponential variation of Tk with Ed one can 
see a semiquantitative agreement between both results 
in Table 1. The remaining quantitative discrepancy can 
be ascribed to the difference between T^ and A as a 
measure of the Kondo temperature. In fact from renor- 
malized perturbation theory one obtains T^/A = 0.746 
in the extreme Kondo limit while this ratio increases be- 
yond 1 when valence fluctuations are allowed^ and the 
Wilson ratio decreases i2£r— 

The NCA has the drawback that it does not fulfill 
Fermi liquid relationships. For example, the spectral den- 
sity at the Fermi energy pd(0) at temperatures well be- 
low Tk differs by about 10 or 20% in the Kondo regime 
from the value predicted by the Friedel sum rule^ A 
detailed comparison of pd(w) between NRG and NCA in 
the Kondo regime is given by Fig. 10 of Ref. HH. One 
can see that in addition to the larger value of pd(0), the 
spectral density is more asymmetric for the NCA. This is 
probably the main reason of the factor near five between 
the magnitude of the dip in S(T) near Tk calculated 
with NCA with respect to the NRG result. Part of the 
discrepancy is probably also due to lack of resolution at 
finite energies within the NRG. For example in models 
of two quantum dots, the split Kondo peaks in the spec- 
tral density are considerably broadened, losing intensity 
(Fig. 11 of Ref. Ull). Another example is the plateau 
in the conductance G(T) observed at intermediate tem- 
peratures T in transport through Ceo molecules for gate 
voltages for which triplet states are important ) 52 ! 53 which 
was missed in early NRG studies, but captured by the 
NCAj 54 ' 55 More recent NRG calculations using tricks to 
improve the resolution^ have confirmed this plateaui^ 3 - 

In any case, the one-level SU(2) limit is the worst case 
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FIG. 2. (Color online) Total occupancy as a function of Ed 
for S —> +oo and different approximations. 




FIG. 3. (Color online) Thermopower as a function of temper- 
ature for E d = — 4F and different values of the level splitting. 



for the NCA, because the real spectral density tends to 
be symmetric, while the NCA improves with increasing 
N for SU(N) models^ 

The absolute value of S is exaggerated within NCA 
for T — > 0. However, as shown in Fig. there is a 
good agreement between the occupancies calculated with 
NRG and NCA. This suggests that using the occupancies 
calculated with NCA and Eq. (|12p a semiquantitatively 
correct result for the linear part of S(T) as T —> can 
be obtained. Instead, while the SBMFA satisfies Fermi 
liquid relationships, the occupancies are not accurate for 
5 — > +oo. Even perturbation theory up to second order 
in V v neglecting spin flip leads to a better result for n, 
although (in contrast to SBMFA) this approach is unable 
to predict the right magnitude of Tk ■ 



B. Temperature dependence in the general case 

In Fig. [3] we represent the NCA results for the Seebeck 
coefficient S(T) as a function of temperature for differ- 



ent values of 5, which represents the splitting due to a 
Zeeman term for example. The spectral density Pd(y) 
for finite 6 has been studied before. 24 i 28 At 6 = 0, /9d(w) 
has a peak just above the Fermi energy with width of 
the order of 2T K , where T K is the Kondo temperature 
of the SU(4) limit (5 = 0. For the parameters of Fig. 
T K is of the order of O.OlXi^ As 5 increases above T K , 
the Kondo peak splits and another peak at an energy of 
the order of 6 above tp appears. This peak above tp 
originates a dip in S(T) for T ~ 6. Note that peaks in 
Pd(w) above ep lead to a positive contribution to 1\ [see 
Eq. l[5])]. which in turn lead to a negative contribution 
to S(T) at the corresponding temperature [see Eq. ([2])] 
Most of the spectral weight of the spectral density lies in 
the charge-transfer peak at energy Ed, which lies below 
the Fermi energy. Thus, when the temperature reaches 
values of the order of the charge-transfer energy ep — Edj 
the thermopower becomes positive [Iq > 0, I± < in Eq. 
©]• 

As discussed above, for 6 — > +oo (the case of one SU(2) 
doublet), S(T) shows one dip at Tk ~ 10~ 6 r and a peak 
near the charge-transfer energy. For S — 0, these qualita- 
tive features remain, but the Kondo temperature is four 
orders of magnitude larger, and the dip near T K is more 
pronounced, due to the larger asymmetry of the peak 
in the spectral density with respect to £ft^ leading to a 
larger integral I\ [see Eqs. {2]| and {3])]. Since Tk changes 
with 5, as discussed below, the dip at Tk displaces to- 
wards larger temperatures as 8 decreases. Instead, the 
charge-transfer peak remains approximately at the same 
temperature and decreases in magnitude due mainly to 
additional broadening of the corresponding peak in the 
spectral density p(w) as the SU(4) limit is approached, 
and also due to some transfer of the spectral weight to 
the Kondo peak. While this transfer is not large, the 
Kondo peak has a larger weight in the integrals I n due 
to the factor of the derivative of the Fermi function [see 
Eqs. © and ©]. 

For S > T K , for example S = 0.02T in Fig. |3J an addi- 
tional dip develops due to the peak near S in pa(uj). As S 
increases further, the dip moves to higher temperatures, 
as it is apparent for 6 = 0.05 and 0.1 in the figure. The 
relative minimum of S(T) near the dip lies at tempera- 
tures of the order of 6 but smaller, probably because of 
the large intensity of the peak near the charge-transfer 
energy, which pushes this minimum to lower tempera- 
tures. For S — 1 this dip turns to a shoulder at the left 
of the charge-transfer peak and for larger S, the dip and 
the peak cross, interchanging the order of temperatures 
for which they appear. 

In Fig. U we show the thermopower as a function of 
T/Tk{S), where Tk{$) is the Kondo temperature for each 
value of 6. Here we determined Tk{8) from the require- 
ment that the conductance G{Tk{5)) = e 2 /h for sym- 
metric leads (Vl = Vr). To a good degree of accuracy, it 
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FIG. 4. (Color online) Thermopower as a function of T/Tk 
for Ed = — 4r and different values of the level splitting. 



FIG. 5. (Color online) Conductance as a function of T/Tk 
for Ed = — 4r and different values of the level splitting. 



is given by the expression 



T K {6) = { {D + S)Dexp 



TtE 1 

~2A 



5 1 
4 



1/2 



6 

2' 



(13) 



obtained from a simple variational wave function] 28 ! 44 
With increasing 8, Tk stays roughly constant until 8 > 
T K = Tr-(O) and then it decreases strongly. From Fig. 
[4] it is apparent that the dip at smaller temperatures re- 
mains at T ~ Tk{8) for all values of 5. One can also see 
that the magnitude of the dip and (the absolute value of 
the thermopower for T ~ Tk) increases as 5 decreases, 
being maximum at the SU(4) point 8 = 0. Note that 
while the magnitude of this dip is exaggerated by the 
NCA for 8 —> +oo, as explained in the previous sec- 
tion, we believe that this is not the case for 8 = 0, be- 
cause the spectral density is naturally asymmetric in this 
case, and since the NCA is a 1/N expansion, its accu- 
racy improves with N in SU(N) models. 8 Furthermore 
NCA calculations of the thermopower of Ce compounds, 
in which orbital degeneracy is important, compared well 
with experimental^ 

In Fig. [5] we show the conductance G(T) calculated 
with the NCA for different values of 8. G(0) is slightly 
larger for 8 — > +oo and lower for 8 — with respect to 
the correct values due to the failure of the NCA spec- 
tral density to satisfy Friedel sum rule^ In spite of this 
shortcoming, the overall shape and temperature depen- 
dence of G in these limits agree with those obtained using 
NRGi 1 -^ 

The point that we want to stress here is that although 
the finite energy features for T ~ 8 and ep — Ed are 
present not only in S(T) but also in G(T), they are much 
weaker in G(T). Thus, the Seebeck coefficient might be 
the appropriate tool to study the electronic structure of 
the system at intermediate energies. 
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FIG. 6. (Color online) Coefficient of the linear dependence 
of the thermopower as a function of the level splitting, for 
E d = -4T. 



C. Dependence on splitting for T —5- 

In this subsection, we present results for the term linear 
in temperature T of the thermopower, as T — > 0, using 
two techniques: SBMFA and the Fermi liquid expression 
Eq. (fT2j) with occupations calculated with the NCA. In 
Fig. |5] we represent —S/T as a function of 8. In the 
SU(4) limit 8 = 0, both techniques agree and indicate 
a very large absolute value of SA/T. In fact, in the 
extreme Kondo limit of an SU(N) model, one has rij = 
1/N and from Eqs. ©, ©, (HUJ), CP and a Sommerfeld 
expansion one obtains 



S 



7T 2 T 



3A 



^sin(27r/iV), 



(14) 



and |SA/T| reaches its maximum value 7r 2 /3 = 3.29 for 
N=4. The value in Fig. [6] is slightly smaller due to some 
degree of intermediate valence. 

As 8 increases, there is little variation until the Kondo 
temperature of the SU(4) limit T K (~ 0.01T in the figure) 
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FIG. 7. (Color online) Occupation of the two doublets rn 
a function of the level splitting, for Ed = — 4r. 



is reached. For larger S, SA/T falls due to the change 
of occupations (see Fig. [7]): the lower doublet becomes 
more populated, while the occupation of the higher one 
decreases, keeping a total occupation slightly below 1. 
While the trend of the curve is the same for both ap- 
proaches, NCA and SBMFA, the absolute value of the 
thermopower decreases too much within the SBMFA for 
6 > lOT^-. This is due to the fact discussed in Section 



IIV Al that the occupation predicted by the SBMFA of the 
lower lying doublet is too large for 6 — >■ +oo. While this 
shortcoming affects the conductance or thermodynamic 
properties in a few %, the effect of this increase is more 
dramatic in the thermopower. 

Instead, since the NCA occupations agree with NRG in 
the limit S —> +oo, the approach that combines NCA oc- 
cupations and Fermi liquid relationships is reliable in this 
limit. The agreement with SBMFA and general expecta- 
tions for the SU(4) model in the Kondo limit, indicates 
that this approach is also reliable for S = 0. 



V. SUMMARY AND DISCUSSION 

We have calculated the thermopower of a model that 
describes electronic transport through quantum dots in 
C nanotubes, in which the disorder does not play an es- 
sential role, and Si fin-type field effect transistors in the 
Kondo regime. This regime can be easily controlled by 
applying a gate voltage that modifies the energy of the 
localized levels Ed- Without disorder and applied mag- 
netic field, the model has SU(4) symmetry, as explained 
in Section HI In this case, the thermopower as a function 
of temperature S(T) has a dip (negative S) at tempera- 
tures near the Kondo temperature Tjc, and a peak (pos- 
itive S) at temperatures near the charge-transfer energy. 



For SU(N) symmetry, varying N with width of the Kondo 
resonance ~ 2Tk fixed, N=4 is the most favorable case 
to have a large thermopower at temperatures lower that 
the Kondo temperature. 

When the SU(4) symmetry is broken, leading to an en- 
ergy splitting S between two SU(2) doublets, by a simple 
symmetry breaking field (like a magnetic field), a new dip 
appears in S(T) at temperatures of the order of S. While 
this dip the peak for positive S at the charge-transfer 
energy are clearly displayed in S(T), the corresponding 
features in the conductance G(T) at equilibrium are very 
weak. This suggest that the study of the thermopower 
might be a more useful tool to study the electronic struc- 
ture at finite energies. Another alternative is to study 
the conductance out of equilibrium from which peaks at 
a bias voltage of the order of ±<5/e are expectedj 30 ' 44 

As 5 increases, the characteristic energy scale Tk(S), 
which determines among several physical scales, the 
width of the Kondo resonance and the temperature at 
which G(T) falls to half its value at T = 0, decreases 
following a simple expression Eq. (|13p . For T <C Tk(S), 
S is linear in T with a coefficient that stays approxi- 
mately constant for 5 < Trt(O). As 5 increases further, 
STk(S)/T for T — > decreases by an order of magni- 
tude but Tk{$) decreases by nearly three orders of mag- 
nitude, so that the linear coefficient increases. For this 
calculation, an approach that combines Fermi liquid rela- 
tionships with occupation numbers calculated with NCA 
gives more reliable results that the SBMFA. 

Because of the nature of the underlying SU(4) symme- 
try, the couplings of the four states involved to the con- 
ducting leads are proportional, and independent of the 
state (Vl to the left and Vr to the right), This fact which 
simplifies the calculations is no more true for weak mag- 
netic field in C nanotubes with disorder,^ 9 and in general 
in multilevel quantum dots^i In these cases, as well as in 
situations with total or partial destructive interference, 
such as transport through moleculeS ) 43 ' 44 i 58 the present 
formalism does not apply and a full non-equilibrium for- 
malism seems necessary^ 3 - In addition here, the total oc- 
cupation of the dot is below 1, in contrast to models with 
multilevel systems in which even occupation is allowed. 53 
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